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ABSTRACT 

We analyze the exchange of dark matter between halos, subhalos, and their envi- 
ronments in a high-resolution cosmological TV-body simulation of a ACDM cosmology. 
At each analyzed redshift z we divide the dark matter particles into 4 components: 
(i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster 
halos, and (iv) the background outside of virialized halos. We follow the time evolu- 
tion of the mass distribution and flows between these components and provide fitting 
functions for the exchange rates. 

The exchange rates show gradual evolution as z decreases to 2, and become more 
steady thereafter. For z 55 2 about 15% of the isolated galactic halos cluster per Gyr 
to become subhalos and a similar fraction of their mass returns to the unvirialized 
background. Mass accumulation onto subhalos is equally shared between previously 
isolated halos and unvirialized matter, and is dominated by accretion from the host's 
diffuse matter beyond z ~ 1. This accumulation is balanced for z ~ 0.5 by subhalo 
disruption at a rate of about half of their mass per Gyr. The diffuse component 
in host halos is built by accreting isolated halos and un-virialized material in mass 
shares of 40% and 60%, respectively, and at z < 0.5 also by disruption of subhalos. 
The unvirialized IGM is enriched mostly by stripping of isolated halos, and at z < 1 
also by mass loss from groups and clusters. 

We go on to use our derived exchange rates together with a simple recipe for 
metal production to gauge the importance of metal redistribution in the universe due 
solely to gravity-induced interactions. This crude model predicts some trends regarding 
metallicity ratios. The diffuse metallicity in clusters is predicted to be ~ 40% that in 
isolated galaxies (~ 55% of groups) at z = 0, and should be lower only slightly by 
z = 1, consistent with observations. The metallicity of the diffuse media in large 
galaxy halos and poor groups is expected to be lower by about a factor of 5 by z ~ 2, 
in agreement with the observed metallicity of damped Lya systems. The metallicity 
of the background IGM is predicted to be (1 — 3) x 10~ 4 that of z = clusters, 
also consistent with observations. The agreement of predicted and observed trends 
indicates that gravitational interaction alone may play an important role in metal 
enrichment of the intra-cluster and intergalactic media. 

Key words: galaxies: formation - galaxies: evolution - cosmology: theory - cosmol- 
ogy:dark matter 



1 INTRODUCTION 



"isolated" galaxy-size halos, while cluster, group galaxies, 
and even satellites of massive galaxies, can be identified 
with "subhalos" embedded in the background of a larger 
1Q 13 /i _1 M ) "host" halo. 



Under the standard assumption of hierarchical galaxy for- 
mation, galaxies reside within virialized dark-matter (DM) 
halos. In this picture, field galaxies can be identified with 
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W e can make a further association by identifying the lock et al. 2 000). For example, we present elsewhere (Kolatt 



intra-cl aster medium with the diffuse matter of massive et al. 2000) an analysis of collision rates of sub-structure 



"host halos" — the mass between the subhalos, and the low- 
density intergalactic medium (IGM) with matter located 
outside of any (galactic-sized or bigger) halo . 

Several phenomena suggest that matter is continuously 
exchanged between grouped galaxies, field galaxies, and dif- 
ferent regions of the IGM. These incl ude cooling f lows in 



Fabian 



cluster s, groups and elliptical galaxies ( |Fa.bian 199j 
& NuIsmi 1994), and galaxy interactions (Barnes & Hern 



quist 1 )92). This interchange of matter between different 



populations is qualitatively expected within the CDM sce- 
nario, where halos are continuously interacting and merging 
into larger halos. However, the analytic descriptions of this 
process, such as the Press-Schechter formalism (PS) and its 



variantsj Press fc Schechter 1974 ; Lacey fc Cole 199^ ; Lacey 



& Cole 1994), fail to match the observed phenomena in two 
important ways: (1) they do not follow the halos as distinct 
entities once they are incorporated in larger halos, and (2) 
they treat accretion into halos as a one-way process, ignor- 
ing expulsion back into the diffuse media due to heating and 
tidal forces. Evidence for the existence of such phenomena 
and the survival of subhalos is provid ed, for example, by 
cluster galaxies' velocity profiles (e.g., <\mram et al. 1994) 
and halo truncation observed via high-resolution de nsity re- 



const ructmnsby grjivjtati^^ Natara- 



jan et il. 1998; Tyson, Kochanski, & Dell' Antonio 



Another indication for expulsion of matter from subhalos 
comes from the relatively high metallicity observed in the 
hot diffuse matter in clusters. 

One way to improve analytic treatments is by semi- 
analytic modeling, in which complex processes (such as gas 
cooling, star formation, and supernova feedback) in galaxies 
within merging dark matter halos are followed via simpli- 
fied recipes. As argued below, there are several limitations 
to current semi-analytic models. In particular, the nonlin- 
ear substructure can be properly resolved only via full-scale 
cosmological simulations, which also provide the spatial in- 
formation missing in semi-analytic models. 



Th e usual semi-analyt i c approach (e.g ., Kauffmann 
White, "~ 



Guidcrdoni 1993; Cole et al. 1994 



Bomerville & 



Primac c 1999) using extended PS merging trees does not 
take into account mass loss and exchange between inter- 
acting halos. Certain aspects of galaxy formation may be 
severely affected by these missing features. For example, dif- 
ferent accretion rates onto clustered galaxies and field galax- 
ies should affect their relative star formation rates. In addi- 
tion, the transfer of material that has been "processed" in 
galaxies into the diffuse media may act to transfer metals 
and heat into the unvirialized IGM (see below.) 

Only recently has the dynamic range of iV-body simu- 
lations become wide enough to allow resolution of h alo sub- 
structure in a cosmological statistical sample (e.g., Klypin 
et al. 1! 99). By utilizing such iV-body experiments it is pos- 
sible to investigate the properties of the hierarchy of halos 
and diffuse matter ([Kolatt et al. 199S|; [Bullock 199s|; [Bui ' 



Alth mgh the IGM is generally associated with a gas compo- 



In this paper, we study the exchange of matter us- 
ing such a simulation. We quantify the flow of dark mat- 
ter among four components: isolated halos, subhalos, the 
diffuse media of host halos, and the unvirialized back- 
ground. These categories roughly correspond to field galax- 
ies, grouped galaxies, the intra-cluster medium, and the 
background IGM. The aim is to add a solid quantitative 
result to the general expectation of matter exchange in the 
hierarchical picture, and also to identify the important ex- 
change processes which might provide input for future mod- 
eling. 

An important example of a process where our results 
are relevant is the large-scale redistribution of metals in the 
universe caused solely by gravitational interactions. In the 
second half of the paper we demonstrate how our derived 
exchange rates shed light on this issue using a crude but ex- 
plicit model for metal enrichment. Our model assumes that 
supernova winds spread processed galactic gas from the disk 
uniformly throughout galactic halos (galaxy-mass subhalos 
and halos) . We associate a fraction of each unit of dark mass 
with gas and assign to it a metallicity in proportion to the 
time it spends inside a galactic halo. By following the flow of 
this "enriched gas", we estimate the effect of gravitational 
exchange on enriching the IGM and the diffuse media of 
clusters, groups, and massive galaxy halos. We make pre- 
dictions for the relative abundance of these populations and 
study how the metallicities should evolve with redshift. 

In we provide a brief description of the simulation, 
the halo finder, and the construction of the halo hierarchy. 
In §|l| we quantify statistically the matter exchange rates 
between the four components and discuss the origin and 
destiny of the matter in each component. In §^ we address 
the evolution of metallicity in the different components and 
compare to observational measurements. We discuss our re- 
sults and conclude in 95l 



2 SIMULATED HALOS 



nent, we are making the assumption that this diffuse gas is tracing 
the diffuse dark matter. 



We u sed the ART code (Kravtsov, Klypin, & Khokhlov 
1997) to simulate the evolution of collisionless DM in the 
"standard" ACDM model (Cl m = 1 - f2 A = 0.3; H = 
WOh = 70 km s _1 Mpc -1 ; a s = 1.0). This model universe 
has a present age to = 13.5 Gyr. The simulation followed 
the trajectories of 256 3 particles within a cosmological pe- 
riodic box of size L = 60 fe _1 Mpc from redshift z = 40 to 
the present. A basic 512 3 uniform grid was used, and six 
refinement levels were introduced in the regions of highest 
density, implying a dynamic range of ~ 32, 000. The formal 
resolution of the simulation is thus / ros ~ 2h~ kpc, and the 
mass per DM particle is m p ~ 1 x lO 9 ft- 1 M . We analyze 
15 saved outputs at times between 2 = 5 and z = 0. 

The identification of halos is a key feature of the analy- 
sis; we try to make it objective and self-consistent, following 
the evolution involving halo interactions and mergers. Tradi- 
tional halo finders utilize either friends-of-friends algorithms 
or overdensities in spheres or ellipsoids to identify virial- 
ized halos. These algorithms fail to identify sub-structure 
with well-defined attributes and errors. We therefore have 
designed a new hierarchical halo finder, based on the bound 
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density maxima (BDM) algorithm ( Klypin et al. 199S). The 
deta ils of the halo finde r are described elsewhere (Bullock 



199S; Bullock et al. 2000), and we summarize below only its 



main relevant features. 

After finding all density maxima in the simulation, we 
unify overlapping maxima, define a minimum number of 
particles per halo (A™ m ), and iteratively find the center 
of mass of a sphere about each of the remaining maxima. 
We compute the spherical density profile about each cen- 
ter and identify the halo virial radius R V [ T inside which the 
mean overdensity has dropped to a value A v i r , based on 
the spherical infall model. For the family of flat cosmologies 
( Om + Oa _ 1), the val ue of A v j r can be approximated by 

(187r 2 +82a;-392; 2 )/(l + s:), 
Sl(z) — 1. In the ACDM model used in the current 
vir varies from about 180 at z 1 to A v i r ~ 340 at 



( |Bryan fc Norman 1998[ ) A 
where x \ 
paper, A 

z — 0. If an upturn occurs in the density profile inside 7? v i r , 
we define there a truncation radius Rt. 

An important step of our procedure is the fit of the den- 
sity profile out to the radius min(i? v j r , R t ) with a universal 



functio nal form . We adopt the NFW profile (Navarro, Frenk 
& Whh]e 1995|), 



PNFw(r) 



(r/R s )(l + r/R s y 



(1) 



with the two free parameters R s and p s — a characteristic 
scale radius and a characteristic density. This pair of param- 
eters could be equivalently replaced by other pairs, such as 
R s and R v i T . Using this fit, we iteratively remove unbound 
particles from each modeled halo and unify every two halos 
that overlap in their _R S and are gravitationally bound. Fi- 
nally, we look for virialized regions within R s of big halos 
to identify subhalos near the centers of big host halos, e.g., 
mimicking cD galaxies in clusters. The minimum halo mass 
corresponding to iV™ m = 50 particles is ~ 5 x lO lo /i -1 M0. 
At this minimum mass though, the finder is incomplete and 
fails to identify some halos. Completeness (i.e., 100% halo 
identification) is reached at M ~ 1.5 x 10 11 hT 1 Mq (cf. Sigad 
et al. 2000). The modeling of the halos with a given func- 
tional form allows us to assign to them characteristics such 
as a virial mass and radius, and to estimate sensible errors 
for these quantities. 

The classification scheme used in the rest of the paper is 
as follows. A subhalo is a halo whose its center lies within the 
virial radius of a larger halo A host is a halo that contains 
at least one subhalo. An isolated halo is any halo that is 
not a subhalo and is also not a host. Finally, in § ^| the 
combined set of subhalos and isolated halos are identified as 
galactic halos. All of the mass that is not contained within 
any identified halo is referred to as unvirialized. Note that in 
the PS formalism subhalos are not taken to be independent 
quantities, and that without a lower mass cutoff for halos, 
the unvirialized component is not well defined. 



t "Larger" means the host is at least 25% more massive than the 
subhalo. In the current application we limit ourselves to the first 
level of subhalos, but the classification scheme can straightfor- 
wardly be extended to deal with many levels of subhalos within 
subhalos. 



3 MATTER EXCHANGE RATES 

At any given output time, we assign each mass particle to 
one of four components, and label each of them as follows: 

(i) Isolated halos (I) 

(ii) Subhalos (S) 

(iii) Diffuse matter in host halos (D) 

(iv) "Unvirialized" diffuse matter (U) 

In the next section we will relate the first two components to 
galaxies and the last two (diffuse components) to intergalac- 
tic gas, but in this section we simply study the exchange 
rates of mass among the four components. 

The division into components clearly depends on the 
details of the halo finder, and in particular on the minimum 
mass imposed. For example, had this mass been set to be as 
small as the particle mass, all the particles would have been 
associated with "virialized halos." Possible incompleteness 
of our halo finder near the minimum mass may slightly af- 
fect the subhalo population but this is not a major concern 
here because the subhalo mass function is somewhat flatter 



than the distinct halo mass function ( Bigad et al. 200C ) such 
that most of the mass is in subhalos more massive than the 
completeness limit of 10 11 h -1 Mq. 

At every output time, we compute the total mass in each 
component as well as the rates of mass exchange between 
the components, all per unit volume. We denote the average 
density corresponding to component x by p x , and the rate 
of density flow from component x to component y by p xy . 
The net flow between the two components is thus p xy — p yx . 
The exchange rate is estimated by counting the amount of 
mass gained and lost between each component and dividing 
by the elapsed proper time between outputs. The timesteps 
are selected to be typically separated by ~1 Gyr. In the 
following, the mean densities are referred to in comoving 
units of 10" 2 (p), where (p) = 8.3 x 10 10 /i~ 1 M o ( /i^Mpc)" 3 
is the mean comoving density of the simulated cosmological 
model. The flows are referred to correspondingly in units of 
10- 2 <p)Gyr- 1 . 

Figure ^ is a schematic diagram showing the distribu- 
tion of mass among the 4 components and the flows between 
them, at two different times, z = and 3. The area associ- 
ated with each component is proportional to the mass in that 
component, and the thickness of the arrows is a monotonic 
function of the flow. Figure ^ depicts the redshift evolution 
of the fractional density in each one of the components. 

In presenting the exchange rates we can either focus on 
the origin of the incoming mass to each component or on 
the destiny of the outgoing mass from each component. The 
information content is redundant, because the x origin of y 
is the same as the y destiny of x in the previous time-step, 
but they allow two different angles of view. Figures ^| and 
^| summarize the exchange rates in these two ways. They 
show the evolution of total mass density in each of the four 
components, and the flow rate into (from) each component 
from (into) each of the other components. We note that the 
general variation in time of many of the rates is slow. We 
fit this weak time dependence of the measured rates with a 
quadratic function in z = log(l + z), 



log p = log(po) + bz + az 



(2) 
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Figure 1. Mass distribution among the four components and exchange rates. U=unvirialized background, I=isolated halos, D=diffuse 
in host halos, S=subhalos. Mean densities are quoted in percentiles (i.e., units of 10 — 2 {p)), and flows in 10 — 2 (p) Gyr - 1 . The areas are 
proportional to the mass in the corresponding component. The thickness of the arrows is monotonic with the flow. Left: z = 3. Right: 

2 = 0. 
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Figure 2. Fractional density of each one of the four components 
U,I,D,S (cf. Fig. F|) as function of redshift. 



Table 1 



Density fit coefficients (logpo),b,a* 



I 



D 



U 



10.132 10.154 9.591 10.439 
0.906 0.656 0.497 0.092 
-2.287 -4.218 -3.220 -0.013 

z = log(l + z) 

* log(px) = log(po) + bz + az 2 
\p]=br 1 M e) (b- 1 M V c)- a 



We assign equal weights to the measured rates in the 
different time steps; this is under the assumption that the 



errors are dominated by cosmic variance, which we are not 
trying to model in detail. Table 1 displays the parameters 
for the quadratic functional fits to the density evolution of 
each component, p(z), the fits are accurate to 10 — 15%, 
and Table 2 displays the parameters for the functional fits 
to the whole matrix of flows, p xy . For example, from I to 
D, log p ID (z) = 9.306 + 0.0435 - 1.242z 2 , which equals 8.88 

" at z = and 7.6 x 10 8 at 
3 Gyr -1 . In general, the 
evolution slows down with time because of the ACDM cos- 
mology; with flm = 0.3, the characteristic epoch for loitering 
is z ~ 0.7. 

Table 2 



at 2 = 3; thus p ID = 2.02 x 10 9 
z = 3, in units of ft -1 M@(h~ Mpc) 





Rate fit coefficients (logpo) 


,b,a* 


To 


I 


D 


S 


U 


From 








9.306 


9.027 


9.347 


I 




0.043 


-0.872 


-0.615 






-1.242 


0.694 


1.225 




9.367 




9.305 


9.408 


D 


-1.621 




-1.673 


-0.621 




0.776 




-0.349 


-1.212 




8.961 


9.580 




8.682 


S 


-1.589 


-2.965 




-0.462 




1.436 


0.957 




-0.258 




9.390 


9.492 


8.925 




u 


0.109 


0.233 


-0.842 






1.086 


-1.277 


0.830 





z = log(l + z) 

* \og(p xy ) = log(po) + bz + az 2 
[p]=h- 1 M (h- 1 Mpc)- 3 Gyr- 1 



3.1 Total mass distributions 

The main features of the total mass evolution in each com- 
ponent in the redshift interval z = 4.5 to are as follows: 
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Figure 3. Time evolution of the mass density in each of 
the 4 components (thick black curve in each panel) in units of 
h~ 1 M ( r,( h~ 1 Mpc)~ 3 (right axis) or equivalently 10 -2 X (p), the 
comoving average density of the simulated cosmology (left axis, 
dropping the "Gyr -1 "). Also plotted are the incoming flow rate 
into each component from each of the other components ("Ori- 
gin") in units of 10 -2 X (p)Gyr — 1 (left axis). The symbols are 
the measured flow in the given time interval, and the curves are 
quadratic fits. Isolated halos I=filled circles, dotted lines, ma- 
genta; diffuse matter in host halos D=filled squares, short-dashed 
lines, blue; subhalos S=empty stars, long-dashed lines, orange; 
unvirialized matter U=open squares, solid lines, red. 



1.5 2 2.5 3 3.5 4 4.5 




log(l I z) 

Figure 4. Destiny: same as Fig. but showing the outgoing flow 
rate from each component into each of the other components. 



(i) The isolated component, I, grows from 5 to 19, with 
most of the relative growth occurring at z > 2. 

(ii) The diffuse component, D, grows continuously from 

0. 2 to 19. The growth continues to be significant down to z ~ 

1, long past the time when I component growth stagnates. 

(iii) The subhalo component S grows at roughly a con- 
stant rate until z ~ 1, from 0.2 to 3.6, and then becomes 
rather constant. 

(iv) The background U is depleted slowly, from a density 
of 87 before z = 3 to 53 at z — 0. 



3.2 Origin and Destiny 

The main characteristics of the origin of each component 
are as follows: 

(i) The origin of I (isolated halos) is dominated by ac- 
cretion from U (unvirialized component) at a gradually de- 
creasing rate, from 12 to 2.6 (recall that the flow units are 
10~ 2 (p)Gyr _1 ). About 100% of I is being added every Gyr 
at z ~ 3, while this fraction drops to ~ 20% at * = 0. (Note 
also that ~10% of I comes from D and S; this is mostly a 
consequence of small host halos turning into simple halos 
because of the disruption of all their subhalos). 

(ii) The main source of D (diffuse matter in host halos) is 
accretion from U, at a slowly varying rate, from 1.1 to 3.6. In 
parallel, there is an important supply to D from infalling I, 
which is always about 1/3 of the total incoming flow. About 
100% of D is being added every Gyr at z ~ 2.5, dropping 
to 40% at z = 0. The expulsion rate from S into D, corre- 
sponding to the disruption of subhalos, grows continuously 
after z ~ 3, and becomes higher than the infall rate from I 
into D at z < 0.4, and comparable to the accretion from U 
at * = 0. 

(iii) The origin of S at z > 1 is divided about equally 
between I and U, at roughly a constant rate of ~ 1 each. 
At z < 1, the origin of S becomes dominated by accretion 
from D, reaching a rate of 2.5 at z = 0. The input to S is 
about 100% of S per Gyr at z ~ 2.5, and 70% after 2 = 2. 
It is likely that fraction of the accretion to S from U is due 
to halos that formed outside of a host and subsequently fell 
in, however the time resolution of the simulation outputs 
analyzed did not allow the identification of the intermediate 
I component. The same applies to the fraction of matter 
that went to S through a D phase. 

(iv) The little input to U comes mostly from expulsion of 
I, and at z < 0.5 also from D, at the level of a few percent 
of U per Gyr. 

The main features of the destiny of each component are 
as follows: 

(i) Isolated halos, I expel mass mostly into U, at roughly 
a constant rate of 2.5-3.5. I also turns into D at a growing 
rate that becomes comparable to 2.5 at z < 1. The latter is 
likely due to I halos falling into groups and clusters and then 
being disrupted into D within a single timestep, without 
ever being identified as S subhalos. At z ~ 4.5, about 100% 
of I is being lost (to U) per Gyr, while after z = 2 the 
outgoing rate (to U and D) is about a 1/3 of I. This is as 
expected from the relative masse s of the isolated halos at the 
two redshifts ( Sigad et al. 2000 ), less concentrated density 
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profiles at higher redshift (|Bullock et al. 200C ) , and a higher 
merger rate before z ~ 2 ( |Kolatt et al, 1!)99| ). 

(ii) The diffuse component D outputs mass in roughly 
equal parts to U and I (the latter being mostly D+S turning 
into I). At z < 1, a significant part of the outgoing mass goes 
also to S. The output rate from D to U grows from 0.3 to 3. 
The total output rate is about 250% of D per Gyr at z ~ 4.5, 
and about 1/3 of D after z = 2. 

(iii) The disruption of subhalos causes S to lose roughly 
a constant ~40% of its total mass every Gyr to D, with a 
rate growing from 0.1 to 3. There is also some outflow from 
S to U at roughly a constant rate of 0.36-0.6, and transition 
from S to I at a rate 0.7-1.0. 

(iv) The main destiny of U is I. The accretion from U to 
D becomes comparable after z = 1, at a total accretion rate 
into halos of 6, which is ~ 10% of U per Gyr. 

The largest absolute flow at z > 2 is accretion from U 
to I, at a rate of ~ 10, partly balanced by a flow back from I 
to U at a rate of ~ 4. They reach a near balance by z = 0, at 
a rate of ~ 2.5. At z < 1, the dominant flow is the accretion 
from U to D, at a rate of ~ 4, which is approaching a balance 
with the back-flow from D to U only near z — 0. 

At z > 2, all the halo components (especially I and 
D) absorb mass at a relative rate comparable to their own 
mass every Gyr, and lose mass at a lower rate. At z < 2 the 
fractional mass inflow is typically 10 to 50% of each halo 
component per Gyr, and the outflow rate is only slightly 
lower. The S component exchanges mass at a high rate of 
40-70% of its own mass every Gyr. 

3.3 Net mass exchange rates 

The evolution of the net mass transfer rate between every 
two components (p xy — p yx ) is shown in Fig. |B| Of course, 
the net transfer can be very different from the actual flows 
in each of the opposite directions. For example, when the 
flows in the opposite directions are equal, the net transfer 
is zero even when the entire population of each component 
has been exchanged. The main lessons from the net transfer 
rates are as follows: 

(i) The net exchange U-I at high z is dominated by the 
accretion from U to I; it slows down in time while the expul- 
sion from I to U increases. At z ~ 0.3 the backfiow becomes 
comparable to the build-up of I from U. 

(ii) The net flows from U to S and from I to S are pretty 
constant at ~ 0.25 — 0.35. 

(iii) There is always a significant net flow from U to D, at 
the level of ~ 1. It exceeds the net flow from U to I at z < 1. 
At 2 < 0.5, the evaporation of D to U competes effectively 
with the slowly increasing accretion from U to D, resulting 
in the decreasing net flow. 

(iv) The net flow from I to D starts low at early times due 
to the near absence of D material at z > 3.5. It increases 
gradually until z ~ 0.5 by the clustering of I accompanied 
by tidal stripping. At lower z the signal is dominated by 
artificial fluctuations. 

(v) The net exchange between D and S is very low at 
z > 1, where the accretion from D to S is only slightly higher 
than the tidal stripping from S to D. At later epochs, the 
tidal stripping becomes dominant causing the net exchange 
to reverse its sign and grow to ~ 1 at z = 0. 




Figure 5. Net exchange rate between every two levels (x,y), i.e., 
X — Y = p xy — pyx, as calculated from the fits of Table 2. Solid 
segments are for positive exchange rate while dashed segments 
represent redshift zones where p xy < p yx - 



In summary, we have presented the mass exchange rates 
between several components of cosmological matter in our 
simulations. We confirm the generally expected trend that 
the mass fraction of "unvirialized" background material (i.e., 
all mass outside of any identified halo) falls steadily as a 
function of time, directly or indirectly fueling the mass ac- 
cumulation in isolated halos, subhalos, and the diffuse media 
of host halos. Although the flow of mass is dominated by the 
accumulation of unvirialized background material into and 
onto isolated halos and host halos, we find that there is a 
significant amount of mass loss from halos back to the un- 
virialized component and to the diffuse component of large 
host halos. In the next section we explore one possible im- 
plication of our derived mass exchange rates. 



4 INTEGRATED HISTORY: METALLICITY 
4.1 Model Predictions 

Interesting astrophysical implications may be extracted 
from the integrated history of the mass in the different com- 
ponents. As an example, we describe an attempt to learn 
about the possible role of gravitational effects in determin- 
ing the metallicity of the gas in the different components. 
We make the crude assumption that the baryons trace the 
mass distribution everywhere and at all times, and that the 
star formation rate per unit mass and the resulting metal- 
licity yield are the same in all galactic halos at all times. 
This is based on a scenario in which supernova-driven winds 
are sufficiently energetic to drive the gas out of the galac- 
tic disks and distri bute it in quasi-static equ ilibrium in the 
galactic halos (e.g., Mac Low fc Ferrara 1999 ), but that this 
feedback is not strong enough to drive the processed mate- 



rial out of galactic h alos (e.g., ^ader 1986; 



Ferrara, Pettini 



fc Shechekinov 200C). The hypothesis we try to evaluate is 
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that gravitational interaction is sufficient to strip the en- 
riched matter from galactic halos and distribute i t in the 



Gnedin 



diff use c omponents (David, Forman, & Jones 1991 
199f~] 

This is of course an extreme hypothesis. Because star 
formation rates grow with the baryon density, the baryons 
that become stars are highly concentrated due to gas cooling 
in the centers of halos; and supernova-driven winds may not 
be able to distribute the metals produced by those stars 
throughout the halos. Doubtless in the real universe there 
is a mixture of the purely gravitational processes considered 
here, on the one hand, and the effects of baryonic physics 
(gas cooling, star formation, supernovae), on the other. In 
order to estimate the relative importance of gravitational 
vs. baryonic processes, however, it is simplest to consider 
extreme cases first. 

Based on these simplifying assumptions, we virtually 
assign a fixed fraction of gas to each mass particle of the 
simulation, and associate with it a metallicity which grows 
in time in proportion to the time the particle has spent in a 
galactic halo. As defined above, the galactic halos are all of 
the subhalos plus all of the isolated halos. 

We focus first on the diffuse component in host 
halos (D). We identify host halos more massive than 
10 13 - 6 /i _1 M a = 4 x 10 13 /i _1 M Q with "clusters", and those 
in the range 10 12 - 10 13 - 6 /i _1 M Q with "groups". Although 
we use the term "groups" for this second component, these 
objects, especially those in the lower half of the mass range, 
are probably better associated with massive galaxies hosting 
satellites than the halos of what are typically referred to as 
galaxy groups. This should be kept in mind when comparing 
our predictions to observations, as discussed in §4.2. 

For each host halo, we keep track of the fraction of the 
diffuse mass that spent less than a given amount of time 
in galactic halos. These fractions are averaged over all the 
host halos of each of the two mass classes. Fig. ^| shows the 
average time distributions at two different epochs. At z = 
0.7, when the ACDM universe considered here had an age 
of ~ 7 Gyr, the time distributions for the two mass classes 
are very similar; about half of the diffuse component has 
not spent any time in galaxies, namely it has been accreted 
directly from the unvirialized background. About 5% has 
spent more than 4.5 Gyr (similar to the solar age) inside 
galaxies. At z = 0, the distributions become quite different. 
First, while for the massive hosts the fraction of D material 
that has not spent any time in galaxies is still ~ 50%, this 
fraction is only ~ 20% for the less massive hosts. Second, 
while the D faction of massive hosts that has spent more 
than 4.5 Gyr inside galaxies is ~ 10%, the corresponding 
fraction of less massive hosts is ~ 60%. Thus, the diffuse 
component in "clusters" has not been enriched significantly 
by processed material since z — 0.7, while in "groups" it has 
been enriched significantly during this recent epoch. This 
difference hints at a different metallicity enrichment history 
in clusters versus groups. 

Figure ^ shows the overall metallicity production rate 
in our scheme as a function of redshift. This is estimated as 
the time derivative of the average metallicity, mass weighted 
in all halos. The general behavior of the metal production 
rate is a rise of about an order of magnitude between z ~ 4 
to 2 ~ 1 — 1.5, followed by a gentler decline of about a factor 
of 3 to the rate value at z = 0. The interpretation of such 




log[M/(h-'M s )] 
13.6-15.2 
- 12.0-13.6 



2 4 6 8 10 12 

t (Gyr) 

Figure 6. The distribution of time spent by the diffuse compo- 
nents in "galaxies", at z = (thin) and z = 0.7 (thick). The host 
halos are divided by mass into "clusters" and "groups" . 




Figure 7. Metallicity production rate versus redshift. The sym- 
bols represent time derivatives of the halo mass weighted metallic- 
ities while the solid line is a parabolic fit to the equally weighted 
points (relative units). 



a trend in light of observations will be further discussed in 
the last section (§[|). 

Next, we consider the integrated metallicity in the dif- 
ferent components defined above: the galactic halos I and 
S, the diffuse component D divided by host-halo mass into 
"clusters" and "groups" , and U, the unvirialized IGM. Since 
the absolute values of the yield and the gas fraction are un- 
known, we focus on the measurement of relative abundances, 
between different epochs or different environments. Figure 
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^| shows the computed evolution of metallicity Z x in each of 
these components, in units of the metallicity in clusters at 
2 = 0, Zd(0). 

The average metallicity values of the different virial- 
ized components grow roughly in proportion to time for 
2 ~ 4. In more detail, they all grow between 2 = 4 and 
2 at a similar rate of about AZ/Z d (0) = 0.08 Gyr" 1 . After 
2 ~ 2 (namely, during the last 10 Gyr), the growth rate in 
clusters continues at a similar pace, but the growth rate in 
the other components speeds up somewhat, to AZ/Z c \(0) = 
0.17, 0.20, 0.24Gyr~ 1 for the diffuse component in groups, 
field ga laxies (I), and clustered galaxies (S) respectively. 



Ths average metallicity in the diffuse component in 



clusters and groups (D) is about one half that in the galaxies 
(I and S). This is in agreement with our finding in § ^ that 
the (enriched) flow from I and S to D is roughly comparable 
to the (fresh) flow from U to D. 

The metallicity in clustered galaxies (S) is higher than 
that of field galaxies (I), by about one third. This is because 
most S halos are old; they formed early as I halos, then fell 
into groups and clusters, and thus typically had a long time 
to produce metals. Many of the current I halos are relatively 
young, and therefore less metal rich. 

The faster growth rate in groups versus clusters at z < 2 
leads to a present metallicity in groups almost twice the 
met allicity in clusters. This is in general agreement with the 



diff eren ce in the time distribution seen in Fig. H. 

The average metallicity in the unvirializea (U) back- 
ground is roughly constant since z ~ 4, at a level of 



Zi GM /Z cl (0) ~ 10" 



10" 



4.2 Oomparioon with obacrvationa 



In order to compare our predictions with observations we 
have to consider in some more detail the association of the 
simulated halo components with observed objects. The asso- 
ciation of the most massive host halos with galaxy clusters is 
natural at all redshifts. At low redshift, the low mass r ange of 
host halos fits well the mass range of galaxy groups (Hwang 



et al. 19991 pavis, Mulchaey, fc Mushotzky 1999| ) or even 
massive galaxies with satellite companions. It is less obvious 
what observations to associate with the low mass host halos 
at high redshift. The choice might be damped Lya systems 
(DLAS), or even Lyman- limit systems and high column den- 
sity (~ 10 17 cm -2 ) Lyman-forest clouds. At high redshift, 



shown in the figure); this improves the agreement with the 
observations, for which massive clusters were preferentially 
selected, especially at high redshift. 

The observational estimates of the 2 = metallicity in 
groups and massive galaxies with satellites are less secure, 
and one cannot yet confirm or refute the predicted trend of 
Z gr (0) ~ l.&Z c i(0) ~ 0.5Z Q . The average metallicity of the 
diffuse component in ric h galaxy groups is measured locally 
to be Zgr(0) ~ 0.35^0 ( |Hwang et al. 1999[ ) with the scatter 
being dominated by an intrinsic variance between groups of 
at least 0.1, partly due to Poisson noise and the small num- 
ber of galaxies per group. Davis, Mulchaey, & Mushotzky 
(1999[ ) studied 17 poor groups and found metallicities of 
~ 0.1Z e for 0.7-1 keV groups (~5-8x lO 13 h -1 M ) and 
~ 0.2 — 0.5Zq for more massive groups. These values are 
smaller than our results. 



However recently Buote (2000) showed that for 12 
bright groups, much higher metallicities are obtained if two- 
temperature models are applied with or without cooling 
flows. The group metallicities increased from Z sr = 0.29Z© 
using a single temperature model to 0.65—0.75^0 depending 
on the details of the two temperature models. These values 
are much more consistent with the values presented here. 
Perhaps, though, a better comparison with our small-mass 
hosts is the metallicity estimates of hot gas in massive el- 
liptical galaxies at z = 0. The meas ured values (Matsumoto 



et al. 1997; Buote & Fabian 1998) range all the way from 



Z = (0.19 ± 0.12)Zq (when single-temperature models are 
invoked) through Z = (O.6±O.5)i?0 (when better fits of two- 
temperature models are considered), and up to Z = O.QZq 



for the galaxies of best signal-to-noise ratio ( Buote & Fabian 



199S ) . The estimate for the metal licity of the Galactic halo 



interstellar gas, however is lower (Savage & Sembach 1996) 



and stands on about 0.1 solar for most elements. The ob- 
served metallicities would be consistent with our predictions 
for Z SI (0) / Z c \(0) if a substantial fraction of the less massive 
host halos which we classified as "groups" actually corre- 
spond to the massive galaxies of Z ~ 0.6 — 0.9Z Q . 

If DLAS are the counterparts of our "group" halos at 
high redshift, then the relevant metallicity measurements 
( |Lu et al. 1998; ; [Prochaska fc Wolfe 199S| ) give Z OLA ^ 
0.07^0 at 2 J; 1.5. The measuremen t for a sample of DLAS 
with 2 <; 1.5 (IPettini et al. 199SI) is Z ~ 0.LZ Q . These 



data points are shown in Fig. |8j. When compared to to- 
day's cluster metallicity, the observed values correspond to 
the unvirialized IGM can be identified with Lya systems of Z gT (z ~ 1.5)/Z c i(0) ~ 0.2, in surprisingly good agreement 



very low col umn densities, 



simulations (Dave et al. 1998; Lu et al. 1998) to correspond 



whic h are shown by 



to the mean mass density at that time. 

For clusters today the observational estimate is Z c \ (0) = 
(0.29 ± O.Ol)^©, with a relatively small t rue intrinsic vari- 
ance between clusters of order Q.QQZq (Edge & Stewart 
199l| ; [Yamashita 1992j ). We show in Fig. |^ three data points 



with erro r bars based on measurements of Fe relative abun- 
dan ce by Mushotzky fc Locwenstein (1997 ) (at 2 ~ 0.2,0.4), 
and ponahue fc et al. (2000[ ) (at 2 ~~08). Thus, the ob- 



served metallicity in the diffuse component of clusters does 
not seem to vary significantly between 2 = and 0.4, and 
not to change by more than a factor of two out to 2 ~ 0.8. 
This is in agreement with the model predictions. An even 
weaker evolution rate is predicted by our model for the most 
massive clusters in the simulation, M > lO 14 /i -1 M0 (not 



with the model predictions. If, alternatively, we associate 
DLAS at high 2 with groups and very massive galaxies with 
satellites at 2 = (both have similar hydrogen column den- 
sities in the range 10 20 - 10 21 atoms cm 2 ), we obtain a 
ratio Z gI (z ~ 1.5)/Z gr (0) ~ 0.26, also in agreement with 
the model predictions. However, there are certain caveats 
associated with this comparison. For example, unlike the 
measurements in groups at 2 = 0, the observed metallicities 
at high 2 were not evaluated by Fe abundance. Also, these 
abundances refer in large to the cold gas component, and 
are relevant to our model predictions only if the assumption 
of proper mixing between the cold disk and the hot halo 
is valid. Finally, the large scatter in these measurements at 
high 2 weakens any conclusion drawn on the basis of the 
mean values. 

The predicted metallicity of the unvirialized IGM may 
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Figure 8. Metallicity of the different components as function of 
look-back time relative to the cluster metallicity at z = 0. Here 
those halos whose mass lies between 10 12 and 4 X 10 13 /i — 'Mq 
are called "groups," (thin solid) while those with higher mass are 
called "clusters" (thick solid). The three data points for clusters 
are from the observed Fe abundance in clusters from Mushotzky 
& Loewenstein (1997) (z ~ 0.2 and 0.4) and Donahue et al. (1999) 
(z ~ 0.8). The data point for groups (empty hexagon) at z = 
(shifted for clarity) is from Buote (2000) cooling flow model. The 
two data points for DLAS at z = 1 — 2 are from Prochaska & 
Wolfe (1998) (circle, higher z) and Pettini et al. (1999) (square, 
lower z). The IGM (U) component metallicity is not shown. 



be compared to the very low column density Lya clouds. 
Data for abundances of these exist only for z ~ 2.2 — 4 



(Savaglio 1997 



Dave et al. 1998 



Lu et al. 1998| ). These re- 
sults show metallicities in the range 0.3 — 3 x 10 -4 . Since no 
significant systematic evolution in metallicity is observed in 
this redshift range, these values may be indicative of lower 
redshift values as well. The simulation analysis shows lit- 
tle evolution in the metallicity of the U component, and 
fluctuations in the range Zigm = 0.3 — 1 x 10 _3 Z c i(0) ~ 
1 — 3 x W~ 4 Zq. This agrees well with the observed values. 



5 CONCLUSION AND DISCUSSION 

We have learned that the role of gravity does not end when a 
virialized halo is formed. Gravity is responsible for a contin- 
uous and substantial exchange of matter among the different 
components of halos and diffuse media. This exchange may 
affect the formation and evolution of luminous galaxies in- 
side the halos and their feedback into the environment. It 
should therefore be included in the modeling of structure 
formation. We have analyzed this exchange between 4 basic 
components of DM halos and diffuse media, and provided 
quantitative exchange rates to be used in further investi- 
gations. We then showed, in particular, that gravitational 
interaction and stripping may provide the mechanism which 
produces a significant part of the relatively high metallici- 



ties observed in the diffuse hot gas of big galaxies, groups 
and clusters as well as the unvirialized background IGM. 

A summary of the matter division and exchange rates 
is as follows. Most of the mass within identified halos (?Z 
10 11 h -1 Mq) is located in "isolated" galactic halos before 
z ~ 2. Then about 15% of their mass accumulates in groups 
and clusters per Gyr, while a similar fraction is expelled 
back to the unvirialized background. Subhalos tend to form 
later, mostly by a constant-rate accretion from outside the 
clusters, but at z < 1 they also accrete from the cluster dif- 
fuse medium. The subhalos are constantly disrupted such 
that about one half the subhalo mass is being exchanged 
every Gyr. The diffuse medium in clusters is built by ac- 
creting extra-cluster material, composed on average of 40% 
halos and 60% unvirialized matter. The disruption of sub- 
halos becomes a significant source after z ~ 0.5. Finally, 
the unvirialized IGM is enriched mostly by mass loss from 
isolated galactic halos, and at late epochs, z < 1, also by 
expulsion from clusters. 

The results summarized above highlight important lim- 
i tations of the recipes used by current semi-analytic models 



(|Kauffmann, White, fe Guiderdoni 1993| |Cole et al. 1994 



Somerville & Primack 1999) 
improvement. In particular 



and offer possible routes for 



• Matter expulsion from isolated halos, which is found 
here to be quite significant, is completely ignored in the 
semi-analytic models because they are based on a hierarchi- 
cal clustering formalism which does not take into account 
such a process. Tidal stripping and gravitational heating in 
collisions are known to be more effective in denser environ- 



ments and at higher redshifts (Kolatt et al. 200C). While the 
time dependence can be incorporated in semi-analytic mod- 
els quite straightforwardly, the environment dependence is 
harder to mimic. 

• We notice that "fresh" matter from the unvirialized 
component (U) finds its way to subhalos (S), while current 
semi-analytic models stop the supply of DM and fresh gas to 
halos as soon as they are incorporated in bigger halos, even 
when the subhalo and the host are of comparable masses. 

• We find that accretion from diffuse matter in halos (D) 
onto subhalos (S) becomes the main source of matter for 
S at low redshift, while this process is completely ignored 
in semi-analytic models. Some semi-analytic models do al- 
low for tidal stripping of subhalos, but the matter exchange 
between subhalos and their hosts seems to be a two-way 
process that should be treated as such in the semi-analytic 
models. 

• The fitting formulae Eq. (2) and Table 1 describing the 
exchange rates as a function of time can be used as a first 
attempt to improve semi-analytic treatment of matter ex- 
change. Since we only provide average exchange rates with- 
out specifying their dependence on halo mass, a sensible 
utilization of these formulae is via the implied fractional ex- 
change rates rather than the absolute rates. 

A convenient use of these formulae may be a direct in- 
tegration to yield the density of each component at an ar- 
bitrary time. Integrations of the rate equations in the red- 
shift range of [4 — 0] yield results accurate to about 20% in 
comparison to the actual densities at the upper limit of the 
integration time. The initial conditions were taken to be the 
actual density of the relevant component at the lower limit 
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of the integral (higher redshift). Both the upper and lower 
limits of the integrals can be read directly from Table 1. 

The second part of this paper is a more speculative 
attempt to demonstrate how the matter exchange may af- 
fect measurable quantities. This is an attempt to predict 
metallicities from simulations that do not include explicitly 
any gas, by simply assuming the metallicity of the different 
components is proportional to the time this matter spent in 
galactic halos. We effectively assume that supernova winds 
expel processed matter from the disk out into each halo, so 
that this processed matter spreads uniformly in it. However 
we assume that the SN winds are not powerful enough to 
drive the matter all the way out of individual halos. We ad- 
dress whether gravitational interactions could have acted to 
provide the missing energy and thus to enrich the IGM with 
this material. We avoided any calibration of the specific star 
formation rate and yield in absolute terms, and instead ad- 
dressed relative metallicities, comparing different times or 
different environments. 

We found that the relative metallicities, compared to 
that of galaxies, obtained by gravitational processes in the 
diffuse components of groups, clusters, and the IGM are in 
the same ball park as the observed metallicities at z = 0. 
This scenario seems to predict an anti-correlation between 
the metallicity and the mass of the host halo, which is yet 
to be tested against observations. For example, it may be 
possi ble to measure halo mass es at high z from line pro- 
files flProchaska fc Wolfe 1997|; |Prochaska fc Wolfe 1998j 



groups and clusters. However, close encounters of gala ctic 
halos and the associated harassment ( Moore et al. 199S| ) are 
more fre quent in the crow ded environment of groups and 
clusters (Kolatt et al. 2000) and may thus provide the miss- 
ing ingredient in these models. 

While we mainly addressed the average metallicity in 
the different environments, the scatter about the mean may 
also be of interest. A large sc atter in metallicity at a given 
density has been stressed by |Cen fc Ostriker (1999 ), but 
the mass bins used here can be more easily identified with 
real objects, and thus allow a more direct comparison with 
observations. 

Of course, our attempt to estimate metallicities is ex- 
tremely simplified and it ignores many complexities involv- 
ing gas dynamics and feedback that should be properly ad- 
dressed in more detailed studies. 

For example, the gas cannot be uniform and cannot 
trace perfectly the dark-matter distribution. The cooling 
time is in fact different at different positions in the halo 
and at different times, and it is likely to develop the two 
phases of cold clouds in a hot medium. The feedback may 
not be strong and persistent enough to keep the expelled 
gas hot and transported to large radii in the galactic halos, 
where most of the stripping is expected to occur. 

In reality, the metallicity yield must vary in time and in 
space. A specific worry is associated with the fact that the 
star- fo rmation rate (SFR) is known to drop at late epochs 



z < 



1 dMadau, Pozzetti, fc Dickinson 1997 



Somervillc, Pri- 



Haehncjlt, Stcinmctz, fc Rauch 1997|; |Mallcr et al. 2000D and mack, fc Faber 200C ). This is largely a result of the gas 



thus allow a direct confrontation of prediction and observa- 
tion. The predicted slow evolution of metallicity in clusters 
from z = 1 to is consistent with observations. Match- 
ing the simulated halos with real objects at high redshifts 
is more ambiguous; if DLAS or Lya systems are the high- 
redshift counterparts of local groups and very massive galax- 
ies with satellites, then the predictions agree with the ob- 
served metallicity evolution. Our tentative conclusion is that 
gravitational effects could indeed be largely responsible for 
bringing the metals from galactic halos into the diffuse me- 
dia. 

The process of a gravitational merger mechanism for 
gravitational enrichment has already been addressed by 
Gnedin (1998), in the context of interacting field galaxies 
and the IGM, using < 3 ft _1 Mpc simulations in which sub- 
cell physics was added "by hand". Our N-body results in a 
cosmological volume confirm the feasibility of this general 
idea and demonstrate that gravity may be responsible for 
the redistribution of metals in the diffuse component of clus- 
ters and groups as well. The very low metallicity predicted 
by our analysis for the unvirialized IGM compared to the 
virialized components is in agreement with the predictions 
based on the abovementioned simulations. 

A more elaborate model for the hot gas and metal 
r edistribution in clusters, g r oups, and massive ellipticals 
(|Mathcws fc Brighcnti 1998;; |Brighenti fc Mathews 1999a 



Brighenti & Mathews 19991 ) have demonstrated some of 



the difficulties in combining a mixture of supernovae as the 
source of metals and enrichment of the diffuse media. While 
these models succeed in reproducing the observed abun- 
dances of hot gas in elliptical galaxies, the straight sum of 
the predicted enrichment from the member galaxies fails to 
account for the observed metallicity of the diffuse media in 



processes, such as cooling, that determine the characteristic 
upper limit of galactic masses, which we do not explicitly 
simulate. Furthermore, we do not explicitly take into ac- 
count the effect of gas depletion in old halos on the SFR. 
Two features of our analysis ease this problem. First, our 
analysis does mimic a moderate drop in the SFR at late 
times, if SFR is taken to be the net flow rate (time deriva- 
tive of the density) to each galactic component. This drop 
is due to the combination of three effects: (a) the slow down 
in the growth of fluctuations due to the ACDM cosmology, 
(b) the imposed minimum on the halo masses, and (c) the 
identification of only simple halos with galactic halos. Sec- 
ond, the yield of metals is expected to vary slower than the 
SFR because it results from stars of different ages. 

A related, interesting comparison to observations can 
be made by the examination of the overall metallicity pro- 
duction rate in view of the SFR diagrams as obtained from 
observations (the "Madau diagrams"). A direct comparison 
between the metallicity production rate in our recipe (Fig. 
Q) and recent versions of SFR diagrams (see Somerville et 
al. 2000 for references) indicate two distinct differences: (i) 
The SFR peak is obtained at a higher redshift (z ~ 2.5 — 3) 
than in our recipe for the metal production (z ~ 1 — 2) 
(ii) the drop from the SFR peak value to its z = value 
is slightly bigger (~ 3 — 4) than our value (~ 2 — 3). How- 
ever the link between an SFR diagram and metallicity pro- 
duction diagrams involves the assumption that all metals 
are produced by type II supernovae which closely trace the 
SFR. If in contrast, a large fraction of metals (especially 
Iron, from which the cluster metallicity is calculated) is pro- 
duced in type la supernovae then the one to one mapping 
between SFR and metallicity production becomes less clear. 
Such SNIa contribution to the metal production is now fa- 
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vored by observations (cf. Renzini 2000). In a scenario where 



combined contributions from different SN types dictate the 
metallicity production rate, the observational SFR diagrams 
can be easily reconciled with the metallicity production rate 
as derived by the proposed recipe of this paper. 

Our comparison with metallicity observations is also 
very simplified. For example, our simple model provides the 
volume-average metallicity within the host halos, but in re- 
ality one may expect spatial gradients in metallicity. The 
observations, which are sometimes directed to different po- 
sitions in the halos, e.g., to cluster centers, may therefore 
yield different metallicities. Another complication is that the 
metals may be shared between (at least) two phases of gas, 
cold and hot, as well as the stellar and dust components. 

Our crude treatment of one single "metallicity" is an 
over-simplification, because actual metallicity measurements 
refer to a variety of elements, usually different elements at 
different redshifts depending on the detectability of the cor- 
responding spectral lines. Moreover, since different elements 
in different a bundances are produced through a variety of 



mechanisms ( Mathews fc Brighenti 1998 ; Brighenti & Math- 



ews 19$9a|; prighenti fc Mathews 1999b[ ), the detailed his- 
tory of various processes may be relevant for the metal pro- 
duction. A similar difficulty arises due to the (sometimes 
unknown) ionization state of the medium in which metallic- 
ities are observed, because in our simple model we cannot 
determine the gas temperature. 

Our metallicity estimates should therefore be consid- 
ered very crude and serve only as an encouraging indication 
for the possible important role of gravity in the large-scale 
exchange of metals. 

It will be very interesting to see how important "na- 
ture" (dark matter physics) is compared to "nurture" (bary- 
onic physics). The fact that the Madau-like plot (Figure M), 
based on including only gravitational effects, looks reason- 
ably similar to the observationally based plot is encouraging. 
Semi-analytic modelers should accept the challenge to dis- 
entangle the relative roles of "nature" vs. "nurture." 
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